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SUMMARY 


Several unmanned multiple -target mission opportunities to comets and 
asteroids have been studied. The targets investigated include Grigg-Skjellerup, 
Giacobini-Zinner, Tuttle-Giaeobini-Kresak, Borrelly, Halley, Schaumasse, 
Geographos, Eros, Icarus, and Toro, and the trajectories consist of purely 
ballistic flight, except that powered swingbys and deep space b >rns are employed 
when necessary. 

Optimum solar electric rendezvous trajectories to the comets Giacobini- 
Zinner/85, Borrelly/87, and Tempel(2)/83 and /88 employing the 8.67 kw Sert III 
spacecraft modified for interplanetary flight have also been investigated. 

The problem of optimizing electric propulsion heliocentric trajectories, 
including the effects of geocentric launch asymptote declination on launch vehicle 
performance capability, has been formulated and a solution has been developed 
using variational calculus techniques. Comparison cases have been run on the 
computer and the results have been presented as a paper at an AIAA conference. 

Major improvements have been made to the HILTOP trajectory optimization 
computer program, and a report detailing the changes is being published con- 
currently. 

An er^or analysis of high-thrust maneuvers involving spin-stabilized 
spacecraft has been developed and applied to the Synchronous Meteorological 
Satellite mission. 
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I. INTRODUCTION 


This study of unmanned multiple -target mission opportunities to comets 
and asteroids begins with a presentation of performance requirements for single- 
target comet rendezvous missions via solar electric propulsion. 

Optimum solar electric rendezvous trajectories to the comets Giacobini- 
Zinner/85, Rorrelly/87, and Tempel(2)/83 and /88 employing the 8.67 kw Sert III 
spacecraft r 'ified for interplanetary flight have been investigated and the results 
are presented in Section II and in Appendix A. Launched by either a Titan III E/ 
Centaur or a Shuttle/Transtage combination, the Sert III spacecraft performs 
rendezvous missions to one of the three comets at various points in the comet's 
orbit, after which the spacecraft remains with the comet indefinitely, possibly 
performing scientific measurements throughout the comet's period of revolution. 
Performance requirements and thruster throttling considerations are discussed, 
and an overview of the missions is presented. 

Section III consists of a presentation of several multiple-target mission 
opportunities to comets and asteroids. Each mission investigated involves 
two targets and consists of an all-ballistic trajectory which returns to the 
vicinity of the Earth after having intercepted the initial target. At the first 
passage of the spacecraft by the Earth, a swingby maneuver perturbs the helio- 
centric trajectory either to send the spacecraft directly to the second target 
or to re-target the spacecraft back to the vicinity of the Earth for one or more 
additional swingby maneuvers. The key parameters defining the existence of 
these double-target-via-Earth-swingby solutions are presented in tabular 
form. 

Section IV consists of a solution to the problem of optimizing electric 
propulsion heliocentric trajectories having high geocentric launch asymptote 
declinations. Comparison cases have been run on the computer and the results 
have been published [9] . 
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Section V describes the major improvements which nave been nrade to 
HILTOP trajectory optimization computer program. A report detailing the changes 
is being concurrently published [lO] . The program changes include the incorpora- 
tion of power degradation, housekeeping power, and declination optimization into 
the model. The ballistic swingby-continuation simulation capability has been 
significantly expanded. 

Section VI summarizes the transfer of the QUICKTOP III and CHEBYTOP III 
trajectory optimization computer programs from the NASA Ames Research Center 
to the Goddard Space Flight Center, and their subsequent usage. 

Section VII presents the results of an error analysis of high-thrust 
maneuvers involving spin-stabilized spacecraft, with application to the Synchronous 
Meteorological Satellite mission. In addition to the simulation of the burn, an 
algorithm predicts the expected post-burn errors. The procedure incorporates 
navigational uncertainties in the pre-burn state, and attitude and magnitude errors 
during the burn to estimate errors in the resulting orbital elements. 
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II. SERT III COMET RENDEZVOUS MISSIONS 

Performance requirements for rendezvous missions to three highly 
interesting comets which have perihelia near one astronomical unit from the sun 
have been investigated. Launched by either a Titan III E/Centaur or a Shuttle/ 
Transtage combination, an 8.671 kilowatt Sert III spacecraft, modified for 
interplanetary flight, performs rendezvous missions to one of the three comets 
at various points in the comet's orbit, after which the spacecraft remains with 
the comet indefinitely, possibly performing scientific measurements throughout 
the comet's period of revolution. The comets are Giacobini-Zinner, having a 
perihelion passage on September 5, 1985; Borrelly, having a perihelion passage 
on December 18, 1987; and Tempel II, having perihelion passages on June 1, 

1983 and September 16, 1988. The missions of interest are referenced to those 
perihelia passages. 

The proposed interplanetary version of the Sert III spacecraft was simu- 
lated using the spacecraft model employed by the HILTOP trajectory optimization 

Ci * 1 

program , using parameter values suggested by the Lewis Research Center, 

[ 2 ] 

Cleveland, Ohio . In terms of the engineering parameters used by the HILTOP 
program, the total thrust subsystem, including power conditioners, of the Sert III 
spacecraft has an efficiency rj =. 63376, w'hich includes the 12° thrust canting 
angle of the three operating thrusters; this thrust subsystem accepts a (reference) 
power of 8.671 kilowatts, input to the power conditioners, when all three 
thrusters are operating at full power. The trajectory simulations of interest here 
do not require more than three* full-throttle thrusters operating at one time ; any 
spare thrusters and power conditioners are therefore included in the net space- 
craft mass. The Sert III spacecraft is assumed to have a constant specific im- 
pulse of 2900 seconds; in this approximation, both the efficiency and specific 
Impulse are maintained constant with throttling ratio, whereas in reality they 
vary slightly with throttling ratio. The propellant tankage weight is assumed to 
be ten percent of the propellant weight, and these two weights vary with each 

♦Three is not necessarily the optimum iber. 
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trajectory solution; reserve propellant and tankage must be accounted for in 
the net spacecraft mass. 

The following mass components are assumed to be components of the 
Sert III spacecraft. Three thrusters at 16 pounds per thruster; 48 pounds: 
three power conditioners at 45 pounds per power conditioner; 135 pounds: 
solar array, 340 pounds: solar array structure and gimbals, 20 pounds: 
solar array orientation mechanism, 14 pounds: sun sensor on solar array, 

5 pounds: computer for thruster operations, 10 pounds- and cabling, 7 
pounds. These mass components add up to 579 pounds, or 262.6 kilograms, 
which is the "propulsion system mass" assumed by the HILTOP trajectory 
optimization computer program. All other mass is either propellant, tankage, 
or net spacecraft mass. A switching matrix between power conditioners and 
thrusters is omitted. Variation of the above mass component assumptions does 
not alter the basic results of this study; inert mass component values may be 
juggled between the net spacecraft mass and propulsion system mass, keeping 
the sum constant. 

The thrust vector is assumed to be optimally oriented (along the primer 
vector) at each point of a given trajectory, and thrust angle limits are given in 
Appendix A. The optimization criterion for this study is maximum net space- 
craft mass. Coast phases are optimized. Launch date, departure excess speed, 
and departure asymptote declination are also optimized. The effect of degrada- 
tion of the solar array is ignored. Launches are from the Eastern Test Range, 
with an assumed latitude of 28°. 5; the Titan III E/Centaur is assumed to have 
a range safety constraint corresponding to a parking orbit inclination of 32°. 5, 
and the Shuttle/Transtage is assumed to have one corresponding to an inclination 
of 57°. 

The six standard orbital elements defining the two-body motion of each 
comet In the ecliptic system, in astronomical units and degrees, are given by 
the following table: 
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Glacoblnl-Z Inner Borrelly 

a 

3.516014 

3.610604 

e 

.7075535 

.6242364 

l 

31.87374 

30.32472 

n 

195.20738 

75.26815 

(0 

172.48223 

353.33262 

t 

n 

9/5.6958/85 

12/18.2683/87 


p 


Tempo 1 II 


1983 

1988 

3.035350 

3.036684 

.544894 

.544428 

12.4374 

12.4317 

119.1582 

119.1183 

190.9242 

191.0410 

6/1.5342/83 

9/16.7175/88 


Representative trajectory profiles, projected onto the ecliptic, are dis- 
played in Figs. 1 to 3, for missions rendezvousing at perihelion. The optimum 
launch dates for these trajectories are shown. Missions tend to fall Into multiple- 
year classes, since locally optimum launch dates usually occur each year, when 
the Earth is in Its most favorable position In Its orbit, which may be approximately 
determined intuitively by considering the geometrical configuration of the transfer 
trajectory. All trajectories considered in this study are single-revolution tra- 
jectories, such as those depicted in Figs. 1 to 3, since such trajectories, to 
the comets considered, have minimum solar distances essentially at one 
astronomical unit, which considerably eases the spacecraft thermal design pro- 
blem. All trajectory profiles shown are four year missions, except for one three 
year profile, the Innermost one of Fig. 3. For trajectories which rendezvous 
at the same time but which depart In different years, the longer-fllght-tlme tra- 
jectories will obviously move farther away from the sun than the shorter ones. 

It can be seen that all three comets pass relatively close to the Earth's 
orbtt when the comets are at perihelion; Glacoblnl-Z Inner and Borrelly both 
recede to about six astronomical units from the sun, whereas Tempel II recedes 
to about 4.7 AU. In Figs. 1 and 2, the spacecraft trajectories to Glacoblnl- 
Z Inner and Borrelly have maximum solar distances of about 4 AU, and In 
Fig. 3 the maximum solar distances are approximately In the 3 to 3& AU 
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Figure 2 

COMET BORRELLY RENDEZVOUS MISSION 
USING SERT III TECHNOLOGY 
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Figure 3 

COMET TEMPEL II RENDEZVOUS MISSION 
USING SERT III TECHNOLOGY 
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range. Detailed values of the maximum and minimum spacecraft solar distances 
are given in Appendix A. 

Ftgs. 4 through 6 display the power curves, normalized to the reference 
power, corresponding to the trajectory profiles of Figs. 1 through 3. Gaps In 
the curves represent optimum coast phases. These power profiles represent the 
normalized solar power available to the spacecraft for propulsion along .he given 
trajectory, since this is In accordance with the method by which the HILTOP com- 
puter program simulates trajectories. Housekeeping power Ib not represented, 
and therefore housekeeping power considerations must be included in the net 
spacecraft mass. The HILTOP program trajectory simulation does not employ 
explicit engine throttling, shutdown, and start-up along the powered portion of a 
trajectory, since these operations may be done implicitly, so long as the input 
power to the power conditioners equals the product of the reference power 
(8.671 kw) and the power ratio. Engine throttling considerations are therefore 
accomplished by utilizing normalized power curves such as those shown, and, in 
general, there may be more than one throttling solution available along a given 
trajectory. As a crude example of this, for a spacecraft having three thrusters, 
a power ratio of one-third at a given point along a trajectory may be satisfied, 
in the approximation considered here, by either one thruster at full throttle or 
two thrusters at one -half throttle. The particular solution having the least 
number of thrusters operating at any given instant is usually considered to be the 
best, since this solution implies a higher throttling ratio, which in turn implies 
slightly improved performance in the real-world model. 

If the engines on board the Sert III spacecraft have a minimum throttling 
ratio of one-half, this corresponds to a one-half times one-third equals one- 
sixth, or about .167, normalized power cutoff value, below which the space- 
craft must coast since there is not enough power to operate even one thruster. 
Considering the power profiles corresponding to the four-year missions in Figs. 

4 through 6. it can be seen that this power cutoff value results in imposed 
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coast phases lasting for 800 to 850 days during the portion of flight which Is 
farthest from the sun, and results In an Imposed coast phase of about 300 days 
for the three-year mission of Fig. 6. Due to the HILTOP computer program 
simulation, however, the spacecraft was assumed to thrust during these times, 
the thrust being obviously quite weak; it Is felt that this approximation does not 
greatly affect the results of this performance analysis study, in which the space- 
craft masses, primarily the net mass, are of greatest significance. The space- 
craft masses corresponding to real-world trajectories, having imposed coast 
phases near their aphelia, are therefore assumed to be in good approximation to 
those given in this report, although other mission parameters such as launch 
date may have to be adjusted. Also, some trajectories in this study have very 
brief periods after the launch time, such as in Fig. 5, during which the power 
ratio exceeds unity, a condition which cannot be physically satisfied by the 
assumed 8,671 kw Sert III spacecraft unless a fourth thruster is briefly turned 
on. Such optimum spacecraft trajectories could easily be adjusted, without de- 
grading the performance significantly, so that a maximum of only three thrusters 
would be required. 

The total propulsion system on-time for each trajectory generated is tabu- 
lated in Appendix A; however, these tabulated values merely correspond to the 
duration for which the spacecraft is considered to be thrusting by the HILTOP 
computer program, and therefore are not a realistic assessment of the actual 
thruster on-time, which must be determined by again considering the representa- 
tive power-ratio curves as given by Figs. 4 through 6 together with the 
throttling lower limit of a single thruster. If this lower-limit is assumed to be 
one-half, then, for example, the Borrelly rendezvous power curve of Fig. 5 
may be analyzed approximately as follows: 

Case (1): With Imposed Coast Phase 

In this case, there is an Imposed coast phase from about -1140 to 
about -320 days from comet rendezvous, due to the lack of power available 
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to operate even one thruster. Since there are three thrusters, the pertinent 
values of power ratio are two-thirds — . 667 and one-third - . 333. For 
approximately the first 100 days, 3 thrusters must operate (after which a 
maximum of only two are required for the remainder of the mission), for 
approximately the next 100 days, 2 thrusters are on, and then one thruster 
operates for about 160 days up to the start of the imposed coast phase. Follow- 
ing the Imposed coast phase, one thruster operates for about 150 days and then 
2 thrusters are on for the final 80 days. This amounts to (3 x 100) + (2 x 100) + 
(160) + (150) + (2 x 80) = 970 single-thruster-days, which, if shared equally by 
the three thrusters, implies that each thruster operates for only 323 days. 

Case (2): Without Imposed Coast Phase 

This case is identical to the above except for an additional 820 single- 
thruster -days, which increases the total single -thruster-days to 1790, which, 
when divided by three (thrusters), yields 597 operating days per thruster. 

The Case (1) value is considered to be more representative, but may have 
to be adjusted upward slightly when the imposed coast phase is actually introduced 
into the simulation. The Borrelly-rendezvous power-curve discussed as an 
example is representative of most power profiles for four-year missions involved 
in this comet rendezvous study, and, since thruster on-time values are generally 
less for three year and two year missions, the general conclusion may be drawn 
that the maximum thruster on-time for each thruster of the Sert III spacecraft 
will be "about one year", a situation which will of course improve if a fourth 
"spare" thruster is available to equally share the load. 

Were a fourth thruster to have been assumed in the HILTOP trajectory 
simulation, then the reference power assumed by HILTOP would have been 
4/3 as great as that actually assumed. Consequently, the optimization algorithm, 
sensing a more powerful electric propulsion spacecraft, would have directed the 
launch vehicle to Inject more Initial mass (and, for low power levels, more net 
mass) into heliocentric space at less departure excess speed, requiring the 
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electric propulsion spacecraft to bear more of the burden in effecting the 
rendezvous, and involving some adjustment to the trajectory profile. This pro- 
cess of adding more thrusters, or discrete steps of reference power, Increases 
the net mass up to a point, which would be as close as the discrete reference- 
power variation could get to the true optimum value of reference power, which 
varies with launch vehicle. 

Some of the comet rendezvous mission simulations employed a new theory 
for the optimization of the launch asymptote declination, which is discussed in 
Section IV. Most of the Borrelly missions, some of the Giacobini-Zinner missions, 
but none of the Tempel II missions made use of the new launch optimization model. 
The results of the launch phase optimization are summarized by the tabulations of 
the launch asymptote declination 6 and the parking orbit inclination i in 
Appendix A. 

A performance overview for the Sert III comet rendezvous missions is 
displayed in Fig. 7, in which the top of each bar represents the net space- 
craft mass deliverable to perihelion rendezvous, each notch below the top corres- 
ponds ' » the net mass rendezvous capability at 50, 100, 150, . . . etc. days 
before perihelion, and each bar above the top shows the net mass rendezvous 
capability at 50, 100, 150, . . . etc. days after perihelion. 

The most striking feature of this bar graph is the relatively large net 
masses which the little Sert III spacecraft is capable of rendezvousing with 
Temwel II during the 1988 opportunity. As solar system targets go, these 
ttuee comets are relatively easy to get to, since they pass somewhat near to the 
Earth's orbit. However, as may be seen by examining the comet's orbital ele- 
ments listed earlier in this section, Tempel II is inclined only about 12° to the 
ecliptic, whereas Giacobini-Zinner and Borrelly are in the 30° range. More- 
over, the value of Tempel II's semi-major axis, which is proportional to its 
energy, is closer to the Earth's value than those of the other two comets. These 
facts are considered to be the major contributing factors in the explanation of the 
greater Tempel 11/1988 net mass capability. 
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Figure 7. Performance Overview for Missions to the Three Selected Targets. 
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Considering the TEM/88 TITAN bars of Fig. 7, there is only a small 
penalty In net mass between the 4-year mission and the 3-year mission, 
rendezvousing at perihelion; however, when considering rendezvousing fifty, 
or one hundred, days before perihelion, the penalty is seen to be severe for 
the 3-year mission, but small for the 4-year mission, and this is basically 
because the 3-year mission is more difficult for the little Sert III spacecraft to 
accomplish. Clearly there is a trade-off between mission duration and payload 
delivered, to produce a minimum-cost mission for accomplishing given scientific 
objectives. 

Another obvious fact derivable fi’om Fig. 7 is the lesser interplanetary 
capability of the Shuttle/Transtage relative to the Titan III E/Centaur; what is 
not obvious is the relative launch costs, which hopefully would be less for the 
Shuttle. 

Only the 2-year Tempel 11/1983 mission was investigated, due to the un- 
llkelyhood of a comet rendezvous electric propulsion mission being launched 
before 1981; hence the emphasis on post-perihelion rendezvous, which not only 
yield greater net masses delivered but also are associated with later launch dates. 

Detailed tabular data pertaining to these comet rendezvous missions may 
be found in Appendix A. 
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III. BALLISTIC MULTIPLE-TARGET FLYBY MISSIONS* 


The work described In this section consisted of confirming the existence 
of ballistic multiple -target mission opportunities to comets and asteroids, in 
particular those which pass relatively near to the Earth's orbit, such that one 
or more Earth-swingby maneuvers can be employed to perturb the heliocentric 
trajectory. Details of the Earth swingby technique and complete descriptions 
of the multiple -target mission profiles are discussed by Farquhar, et. al. 

[ll, 12] . 

All missions involve two targets and consist of trajectories which 
return to the vicinity of the Earth after having Intercepted the first target. 

At the first passage of the Earth, a swingby maneuver perturbs the heliocentric 
trajectory either to send the spacecraft directly to the second target or to re-target 
the spacecraft back to the vicinity of the Earth for one or more additional swingbys, 
In which the final swingby sends the spacecraft to the second target. Quite often 
there is more than one way to swing past the Earth at a given Ear-th -pass age, in 
order to produce the desired re-targeting of the spacecraft. 

This work was accomplished using the HILTOP trajectory optimization 
computer program [l, 10 3 . This program generates swingby maneuvers under 
the assumption of the patched-conic approximation, such that the swingby planet's 
sphere-of-lnfluence is assumed to have zero radius as seen from Interplanetary 
space and infinite radius as seen from the planetary vantage point. The passage 
time In the swingby planet's sphere-of-lnfluence is neglected, l.e. , taken to be 
zero In the heliocentric reference frame. 

When, for a particular multiple -target mission, no trajectory solution 
could be found having entirely unpowered swingby maneuvers, powered swingby 
maneuvers were employed, and the minimum AV solution was determined in 
each case by optimizing the post-swingby-leg flight time. Figure 8 depicts the 
relation of unpowered swingby solutions to the wider class of powered swingby 

♦The multiple-target missions described in this section were originally identified 
by R. W. Farquhar of the Goddard Space Flight Center. 
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Figure 8. Relation of Unpowered Swlngbys to Powered Swingbys. 



solutions. The upper curve represents the case In which no unpowered swingbys 
exist, and the lower curve shows how unpowered swingbys occur in pairs, even 
though one or both of the corresponding trajectories may hit the swingby planet. 

Each powered swingby maneuver Is constrained to occur at the mutual 
perlfoct of the approach and departure hyperbolic trajectory segments within the 
swingby planet's sphere of influence; the burn is assumed to be Impulsive with the 
thrust colllnear to the velocity at closest approach. 

The first -target (of a double-target sequence) investigated in this study 
was either the 1977 apparition of the comet Grlgg-Skjellerup (GS/77), the 1982 
apparition of the same (GS/82), or the 1985 apparition of the comet Giacoblnl- 
Zinner (GZ/85). The spacecraft thus launches from the Earth, flies past ore of 
these comets at considerable relative speed, and returns to the vicinity of the 
Earth to perform a swingby maneuver onto a subsequent target. For convenience, 
the following abbreviations are used for three of the comets: 

GS Grigg-Skjellerup 

GZ Giacobinl-Zinner 

TGK Tuttle-Giacobini-Kresak 

An overview of the double-target missions investigated is given in Table III— 1, 
and the six standard orbital elements defining the two-body motion of the targets 
involved, in the ecliptic system, expressed in astronomical units and degrees, 
are given in Table III-2. 

Table III— 3 contains the parameters which characterize the initial helio- 
centric trajectory segments including launch from Earth, flyby of the initial target, 
and arrival back at Earth, for the three cases investigated (GS/77, GS/82, GZ/85). 
In each case, the initial target is approximately crossing the ecliptic at the time of 
encounter, and the trajectory segments are essentially 360° in-ecllptte transfers. 
The trajectory segments defined by Table III-3 apply to all of the associated 
missions to second -targets, summarized by Table III— 1 (except for the unpowered 
swingby to Halley, in which the spacecraft returns to swingby the Earth after three 
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TABLE in - 1 

OVERVIEW OF MISSIONS INVESTIGATED 


First 

Target Second Target 


No. of 
Earth 
Swingbys 


Comments 


GS/77 

GZ/79 

1 

Min AV powered swingby. 


ICARUS/78 

1 

Two unpowered swingby possibilities. 

GS/82 

TGK/84 

1 

Two unpowered swingby possibilities. 


GEOGRA PHOS/ 83 

1 

Three unpowered swingby possibilities, one having relatively 
low flyby speed at Geographos (9.3 km/sec). 


TORO/83 

1 

Possibility of unpowered or min Av powered swingby. 


ICARUS/87 

1 

Min A v powered swingby. 


SCHAUMASSE/84 

1 

Min AV Earth powered swingby followed by Min A V Venus 
powered swingby. 


HALLEY 

1 

Min A V powered swingby (or lone unpowered swingby using 
a three year transfer to Earth). 

GZ/85 

BORR ELLY/87 

2 

First swingby is standoff encounter; second Earth encounter 
has two unpowered swingby possibilities. 


TORO/87 

1 

Powered sw ingby having Earth passage distance 1.2 radii. 


EROS/87 

1 

Family of powered swingbys having low flyby speeds at Eros 
(4 to 6 kni/sec). 


ICARUS/87 

1 

Lone unpowered swingby. 


GEOGRA PHOS/87 

1 

Lone unpowered swingby. 



TABLE III - 2 

TARGET ORBITAL ELEMENTS 


Target 

a 

e 

l 

Cl 

60 

t 

P 

BORRELLY/87 

3.61060400 

.6242364 

30. 32472 

75.26815 

353. 33262 

12/18.26830/87 

EROS/87 

i. 45820548 

. 2227330 

10. 82860 

304. 30050 

178.55830 

5/23.75410/87 

GEOGRAPHOS/83 

1.24472390 

. 3420300 

13.33410 

337.23520 

276.54530 

5/18. 77180/83 

GEOGRAPHOS/87 

1.24495016 

.3420090 

13. 32880 

337.26210 

276.54740 

7/18. 76370/87 

GS/77 

2.96273200 

.6647010 

21.10145 

213.03235 

359.31610 

4/11.00088/77 

GS/&2 

2. 95889800 

.6656651 

21.13392 

213.09221 

359.32122 

5/14. 99629/82 

GZ/79 

3.49071500 

. 7146548 

31.69660 

195.47830 

171. 96939 

2/12. 91951/79 

GZ/85 

3. 51601400 

.7075535 

31.87374 

195.20738 

172.48223 

9/ 5.69583/85 

HALLEY 

17.94348000 

. 9672770 

162. 24046 

58.67151 

111.87092 

2/ 9.39474/86 

ICARUS/78 

1.07785100 

. 8267308 

22.92507 

87.95981 

31.10711 

5/28.72346/78 

ICARUS/87 

1.07795840 

. 8266640 

22.88890 

88.00800 

31.18870 

5/11.87700/87 

SCHAUMA SSE/84 

4.08725300 

.7032880 

11.84050 

80.42600 

57.38400 

12/ 7.36721/84 

TGK/84 

3. 14693400 

.6430245 

9. 93305 

153.73014 

49. 52093 

7/28.41270/84 

TORO/83 

1.36737860 

.4511330 

9.37360 

274.27700 

126. 81560 

12/16. 39710/83 

TORO/87 

1.36736120 

.4511070 

9. 37290 

274.30820 

126.80780 

2/26.47011/87 


Orbital elements obtained from Dr. D. K. Yeomans of Computer Sciences Corporation. 



TABLE m - 3 

PARA METERS CHARACTERIZING INITIAL TRAJECTORY SEGMENTS 


EARTH LAUNCH 

INITIAL TARGET FLYBY 

EARTH ARRIVAL 

Date = 10/30/76 
C 3 = 2. 3 km 2 /sec 2 
6 = -18°. 5 

GS/77 

Date = 4/11/77 
Flyby at 15.2 km/sec 
CD = .20 AU 
CA = 82° 

Date = 10/30/77 

V =1.53 km/sec 

00 

Date « 10/23/81 

2 2 

C =8.0 km /sec 

3 

6 = 9°. 9 

1 

GS/82 

Date = 5/15/82 
Flyby at 15. 3 km/sec 
CD = .37 AU 
CA = 76° 

Date = 10/24/82 

V =2.83km/sec 
00 

Date = 3/10/85 
C 3 = 12.3 km 2 /sec 2 
6 = -4 °. 1 

GZ/85 i 

Date = 9/11/85 
Flyby at 20.6 km/sec 
CD = .46 AU 
CA = 80° 

Date = 3/10/86 

V =3. 51 km/sec 

00 


NOMENCLATURE 


6 -Departure asymptote declination 
CD- Communication distance 
CA -Commun'cation angle (Sun = 0°) 




! 


1 


1 


1 



years instead of one, with the trajectory to GS being the same). Therefore, It 
would be theoretically possible to send several spacecraft, perhaps using a single 
launch, to, say, GS/82 to obtain several different vlewpoltiv. of the comet, and 
each spacecraft could subsequently swing past the Earth differently and continue 
on to a different target. 

Parameters defining the trajectory segments which continue from Earth 
swingby to the final targets are given In Tables HI-4, HI-5, and III— 6, These 
tables, together with Table III— 3, completely define the double-target missions 
which have br an Identified. It Is helpful to refer back to the comments in Table 
III— 1 In order to better understand the latter tables. Multiple swingby p .isslbllttles 
(at the same time) are specified by multiple-values, for certain parameters, 
which are separated by slashes . 

The GZ/85: Eros/87 double-target mission consists of a family of possi- 
bilities having relatively slow flybys of Eros. This family of missions, In which 
the post-swlngby-leg transfer time Is varied, Is summarized by Figure 9. 
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TABLE III - 4 

PARAMETERS CHARACTERIZING FINAL TRAJECTORY SEGMENTS 

FOR GS/77 MISSIONS 


EARTH SWINGBY 


SECOND-TARGET FLYBY 


R =3.21 

AV = 392 meters/sec 

GZ/79 

Date = 2/19/79 
Flyby at 20. 9 km/sec 
CD = 1. 83 AU 
CA = 23° 

R = . .96/31.27 

AV = zero/zero meters/sec 

ICARUS/78 

Dates = (7/15/78)/ (7/24/78) 
Flyby at 29.5/26. 3 km/sec 
CD = . 51/. 65 AU 
CA = 82°/87° 


NOMENCLATURE 

R Passage distance (Earth radii) 

AV Powered swingby incremental velocity 

CD Communication distance 

CA Communication angle (Sun = 0 ) 

(Also applies to Tables III - 5 and III - 6) 
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TABLE III - 5 

PARAMETERS CHARACTERIZING FINAL TRAJECTORY SEGMENTS 


FOR GS/82 MISSIONS 


EARTH SWINGBY 

SECOND-TARGET FLYBY 


R =14.97/1.46 

AV = zero/zero meters/sec 

TGK/84 

Dates = (6/12/84)/ (6/21/84) 
Flyby at 14.5/15.6 km/sec 
CD = 1. 80/1.77 AU 
CA = 43°/42° 


R = 3.07/3.76/12.34 

AV = zero/zero/zero meters/sec 

GEOGRAPHOS/83 

Dates = (3/l6/83)/(7/2/83)/(ll/29/84) 
Flyby at 13. 9/9. 3/12. 5 km/sec 
CD = .10/. 62/1.53 AU 
CA = 120°/66°/39 o 


R =1.38/6.47 

AV = 63/zero meters/sec 

TORO/83 

Dates = (9/29/83)/ (2/13/84) 
Flyby at 10.1/12. 8 km/sec 
CD = 1. 19/. 35 AU 
CA = 65°/90° 


! 

R =1.81 1 

AV = 1014 meters/sec 

ICARUS/87 

Date = 4/1/84 
Flyby at 25.0 km/sec 
CD = 1. 79 AU 
CA = 60° 
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TABLE III - 5 (continued) 


EARTH SWINGBY 


SECOND-TARGET FLYBY 


Earth Swingby to Venus 


R =8.86 

AV = 404 meters/sec 

SCIIAUMASSE/84 


Date = 10/17/84 

Venus Swingby 

Flyby at 12.2 km/sec 

Date = 2/10/84 
R =1.42 Venus radii 

AV = 950 meters/sec 

CD =1.37 AU 
CA = 69° 



HALLEY 

R =17.27/7.10 

Dates = (3/15/85)/ (3/8/85) 

AV = 120/zero meters/sec 

Flyby at 71. 8/78.2 km/sec 
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TABLE III - 6 

PARAMETERS CHARACTERIZING FINAL TRAJECTORY SEGMENTS 

FOR GZ/85 MISSIONS 


EARTH SWINGBY 

SECOND-TARGET FLYBY 




First Passage 


R =64.05 

BORRELLY/87 

AV = zero meters/sec 

Dates = (12/25/87 )/(12/30/87) 

Second Passage 

Flyby at 17. 3/17. 3 km/sec 
CD =.53/.56 AU 

Date = 8/20/87 

R =3.54/3.45 

AV = zero/zero meters/sec 

CA = 125°/122° 



TORO/87 

R =1.20 

AV = 17 meters/sec 

Date = 4/20/87 
Flyby at 15.7 km/sec 
CD = 1. 42 AU 
CA = 44° 



EROS/87 


Date = 5/6/87 

R =5.39 

Flyby at 6 . 3 km/sec 

AV = 220 meters/sec (Min) 

CD = 1. 84 AU 

(see Figure 9) 

CA = 33° 



ICARUS/87 


Date = 6/19/87 

R =13.78 

Flyby at 30.7 km/sec 

AV = zero meters/sec 

CD = .16 AU 


CA = 59° 



GEOGRAPHOS/87 

R =3.79 

AV = zero meters/sec 

Date = 9/26/87 
Flyby at 12. 9 km/sec 
CD = .40 AU 
CA = 90° 
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POST-SWINGBY TRANSFER TIME (days) 
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SWINGBY PASSAGE DISTANCE (ER) 
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IV. LAUNCH ASYMPTOTE DECLINATION OPTIMIZATION 


Introduction. Preliminary performance studies of heliocentric electric 
propulsion missions require some means of correlating initial spacecraft mass 
m o and launch hyperbolic excess velocity V o> In most studies to date, this 
has been accomplished by equating m to the launch vehicle (LV) payload m . , 
which is represented as a non-linear function of the scalar quantity v , the 
hyperbolic excess speed. With few exceptions, this LV payload capability 
assumed has been that corresponding to a due-East launch from the ETR. The 
direction of the launch hyperbolic excess velocity is usually left unspecified and 
is determined as part of the solution to the optimization problem. Using the in- 
direct optimization technique, the solution dictates that V^ be directed along 
the initial primer vector, a requirement that may be in conflict with the 
assumed LV payload capability. 

By properly choosing the point of the (coplanar) injection from a circular 
parking orbit, the geocentric declination 6 of the hyperbolic excess velocity 
may lie within the range 

-i ^ 6 * i, (1)* 

where i is the equatorial inclination of the parking orbit established by the 
launch vehicle. If the launch excess velocity asymptote declination, as deter- 
mined by the solution to the optimization problem, falls within this range and if 
the LV payload capability is compatible with the orbit inclination i, then the 
solution is consistent within the assumptions made and the results are valid. 
However, if 1 6 | > i, then the basic assumptions regaining LV capability are 
in conflict, and it is necessary to formulate the optimization problem to account 
for the dependence of LV payload on direction of the launch excess velocity 
asymptote. 

Although the questionable validity of published high asymptote declination 

solutions has been recognized for some time, no formal treatment of the problem 

♦Equation numbers referenced in this section of the report pertain to this section 
only. 
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has been noted in the literature. The authors had previously developed 11 J a 
technique for adjusting the LV payload to account for the non-coplanar injection 
maneuver required to achieve the geocentric declination of the primer vector, 
which was colinear with V , but this a posteriori correction approach has 
proven unacceptable because the original transversality conditions (TC's) were 
no longer valid. Typically, these TC's resulted in "solutions" which were not 
stationary points. This condition arose because the alignment of V w with the 
initial primer was no longer a necessary condition of optimality, but rather an 
imposed constraint which was in violation of the assumptions used in originally 
formulating the solution. 

In this discussion, a unified treatment of the high asymptote declination 
problem is presented. The LV payload capability is modeled as a function not 
only of the magnitude of V o but also of the inclination of the circular parking 
orbit and of the declination of the launch asymptote. The formulation permits 
the optimization of both the parking orbit inclination and asymptote declination 
or of the asymptote declination subject to a limitation on parking orbit inclina- 
tion to satisfy range safety constraints. The necessary conditions of optimality 
are derived for a typical comet or asteroid rendezvous problem. 

Problem Formulation . High asymptote declinations frequently arise in 
missions to targets that have orbits highly inclined to the ecliptic, such as those 
to certain comets and asteroids. Therefore, we select, for illustrative purposes, 
an optimal rendezvous mission to a single, massless target whose path is defined 
by a specified ephemeris. The extension of the results derived here to other 
missions of interest, such as flybys, orbiters, and multiple-target missions, 
is straightforward. We shall also assume a propulsion system of fixed size u 
terms of mass and reference power Also, overall propulsion system efficiency 
and the specific impulse of the thruster subsystem are assumed given and are 
held constant throughout the mission. 
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The assumed spacecraft and propulsion system models are as described 
in [l] and will not be repeated here* The launch vehicle payload capability is 
assumed to follow the simple exponential law 

m ^ = b i e " (Vc/b2) " b 3’ (2) 


where b , b , and b are pre-determined constants for each launch vehicle 
and V c is a characteristic speed representative of the energy required to achieve 
a specific escape trajectory. For example, for a due- East launch from ETR and 
a coplanar injection maneuver, v c is defined to be the speed inquired at de- 
parture from a low altitude circular reference orbit to achieve a specified hyper- 
bolic excess speed v , i.e., 

00 


v 

c 




( 3 ) 


where v is the circular satellite speed in the reference orbit. Thus, for 
o 

due- East launches and coplanar injection maneuvers, m^ is a function only of 
v^ for a given launch vehicle. Performance data for a large selection of exist- 
ing and potential launch vehicles are presented graphically in Reference [3] 

as a function of v , as defined above, with the reference orbit altitude being 
c 

185 kilometers. The authors have found a least-squares curve fit to the ex- 
ponential law above using, say, 7-10 data points from a given payload curve to 
be a quite adequate and accurate representation of a launch vehicle's performance 
capability. 

To accomodate large launch as; nptote declinations, the same exponential 

law for launch vehicle payload may be used, but the definition of characteristic 

speed must be expanded to reflect the additional energy required to rotate the 

asymptote. This new definition of v is taken to be that given above plus the 

c 

velocity penalties associated with the asymptote rotation. The rotation is 
assumed to be accomplished by first choosing a launch azimuth which establishes 
a given reference orbit inclination i followed by a non-coplanar injection 


33 


i 



maneuver from that circular reference orbit to the desired asymptote, as 
lllusti'ated in the figure on the following page. The velocity penalty incurred 
with non-due-East launches from the ETR is shown graphically in Reference [3] 
as a function of the orbit inclination. This velocity penalty, which we will denote 


Av^ is adequately approximated with a quadratic curve fit of the form 

A W 2+ V + V 


(4) 


Normal range safety limitations restrict the range of inclinations achievable 
through varying the launch azimuth alone. The referenced graph indicates that 
the maximum allowable northerly azimuth will yield an orbital inclination of 
about 48.5 degrees while the maximum allowable southerly azimuth will yield 
an inclination of about 32 degrees. Now, given a reference orbit inclination i, 
it remains to define the velocity penalty Av^ associated with a non-coplanar 
departure from this circular orbit to the desired hyperbolic excess velocity at 
a declination 6 . Assuming the line of nodes of this reference orbit is an open 
variable, one may choose this variable to minimize the angle between the excess 

[4] 

velocity and the orbital plane. This minimum angle is 6 - i. Gunther 
has shown that the minimum incremental velocity required to achieve a given 
v # along an asymptote not lying in the orbital plane from a specified circular 
orbit is obtained from the solution to a quartic equation in the sine of the out-of- 
plane angle. Defining 

s = sin (6 - i); p = v /v ; 

qo O 

2 2 

p=S (p +4); 


2 2 2 
q = S (1 - s ) p ; 


(5) 


- [f( q/2) 2 + (p/3) 3 + q/ 2 ] 3 - [y<q/2) 2 +<p/3) 3 -q/ 2 ] 
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, LAUNCH GEOMETRY 



(A) NON DUE EAST LAUNCH INTO PARKING ORBIT 

(B) OPTIMUM NON COPLANAR BURN OUT OF PARKING ORBIT 





y = 7p 2 /4 - x ; 

W ~ 2 [<°/ 2+ y + ./(P/2-*y) 2 +4(x/2+ /x 2 /4 + s 2 ) ]. 


(5) 

cont. 


then Gunther's solution for the magnitude of the minimum velocity impulse re- 
quired to accomplish the maneuver is 


g 


= V q Jp 2 +3-2 ^(1 + p\v - w 2 )(2 + pw). 


( 6 ) 


and the penalty Av is the difference between v and the velocity increment 
g g 

required if the out-of-plane angle were zero, We. , 


Av = v - ( /v 2 + 2v 2 - v ) . 

g g \V ® 0 C V 


Thus, the definition of the characteristic speed for those cases in which the 
asymptote declination lies outside the interval (1) is 


v = /v 2 + 2v^ + Av. + Av 

c V as O l £ 


= V + V + Av, . 
O g i 


(7) 


Optimality Conditions. The state and adjoint equations for the problem 

under consideration are precisely as formulated in Ref. [l] . The only difference 

in the optimality conditions is the format and content of certain of the trans- 

versality conditions. Specifically, these differences are due solely to the new 

definition of v which is now a function of the direction of V as well as its 
c <® 

magnitude. Whereas before the differential of v was, simply 

d v = (v /v ) d v , 

c ' aa c' 05 

we see from (7) that the equivalent formula now is 

dv = (dv /dv )dv + (dv /d6) d6 
c g * ® g 

+ (dv /S i + d A v./d l) d i 
g i 
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where, from (4) 


dAv { /di = 2 c i i + c 2 , 

and, from (5) and (G) 

dv /d i = - dv /d6 . 
g g 

The derivation of the partial derivatives dv /dv and dv /d6 is straight- 

g 00 g 

forward although somewhat cumbersome. The equations are presented in 
Appendix B. 

After noting that the differential of g, the reference thrust acceleration, 
Is 


dg ~ - (g/m^) dm ( 

where 

dm. - - (b /b ) e '^ v c^ b 2^ dv , 
■C 1 2 c 


and that the differential of V may be written 

00 

dV = (V /v )dv +(n x V )da+ C(V xn)xV / I V xn I]d6, 

where n^ is a unit vector along the North Pole 1 a Is the geocentric right 
ascension of V , one may write the transversal ity conditions associated with 
V by inspection. 

OO 

For optimum launch parking orbit inclination: 

f(dAv /at -dv /d6) = 0. (8) 

‘ g 

For optimum launch excess speed: 

f(dv /dv ) - (A • V )/v = 0 . (9) 

cr oo q oo oo 

For optimum launch asymptote declination: 

f(dv g /d6)-A o - C(V.xn p )xV w /|V w xn p |]-0. (10) 
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For optimum launch asymptote right ascension: 

-A • (ii x V ) - 0 . (11) 

op® 

In the above equations A q is the initial primer vector and f is 

f * [k +k-(l+k)F-gX /m ]dm./dv , 
st t f g -C <■ c 

where k and k are the structural and tan':age factors, respectively, v r is 
s t f 

the final mass ratio, and X^ is the final value of the Lagrange multiplier asso- 
ciated with reference thrust acceleration. 

Equation (11) implies that the right ascension of V b must be equal to, or 

180 degrees from, that of A ;i.e., V must lie in the plane of A and n . 

o 00 op 

If the first term in (10) were zero, which is the result obtained when the effects 
of declination are ignored in the formulation, then one obtains from (10) and 

(11) the familiar result that V must be colineal with A . Usually it is assumed 

° [5] 

that V is aligned with A , however, cases have been found for which the 

oo O 

optimum solution resulted in V being diametrically opposed to A . The fact 

oo O 

that the first te^m in (10) is non-zero means that V w ill be offset from A 

by a finite angle. This offset, as noted above, must be in the plane of A q and 

n^ and, intuitively, we know it must be in the direction of the equator so as to 

reduce v . The amount of the offset of V from A may not be determined 
c ® o 

from (11) as an initial value problem since the variable f is a function of 

variables (v and X ) evaluated at the final time. Thus, 6 must be treated 
g 

as an additional Independent parameter, and (10) becomes another condition to 
be satisfied in the boundary value problem. 

The satisfaction of (8) requires that the term within parentheses vanishes 
since f will normally be a non-zero quantity. Therefore, since the two partial 
derivatives are functions only of initial conditions, one may solve for the l that 
causes the parenthesized term to vanish and thereby eliminate the condition (8) 
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from the boundary value problem. Due to the complexity of the equations defining 

dv /d6, this solution for i must be obtained using an iterative techn'quo. 

S 


The approach to the solution of the problem as formulated above differ^ 
in three basic respects from that of the problem where asymptote declination is 
Ignored: (1) the condition (8) must be solved for the optinv’m parking orbit in- 
clination, given values of v and 6; (2) the asymptote declination 6 must be 

CD 

introduced as an independent parameter and (10) added as an end condition of the 
problem; and (3) the evaluation of V o becomes somewhat more involved. The 
compute tion of V given v , A and 6, proceeds as follows. Denoteas c 

00 00 Q 

the obliquity of the ecliptic such that the matrix 


<J> = 


10 0 
0 cos e -sin € 

0 sin c cos < 


operating on a vector expressed in ecliptic Cartesian coordinates yields the same 
vector in Earth equatorial coordinates. Then the right ascension of the initial 
primer A q may be written 

a. = tan * ["(A cos f - X sin e )/X 1 , 

X i-' yo zo xo J 

where X , X , X are the given ecliptic coordinates of A . Then, the 
xo yo zo ° o 

right ascension of the asymptote is set 


a = or a - a^+n , 


and V is evaluated 


cos acos 6 
sin Ocos 6 
sin 6 

This may be contrasted with the usual definition of V ; 


V = v * J 

® m 


V = v A / j A | . 

Oft Oft o O 
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HILTOP COMPUTER PROGRAM IMPROVEMENTS 


This section describes the modifications and improvements made to 
the HILTOP electric propulsion trajectory optimization computer program [l]. 
New program features include the simulation of power degr; u at ion , house- 
keeping power, launch asymptote declination optimization, and powered and 
unpowered ballistic multiple swingby missions with an optional deep space 
burn. 

Power Degradation . The power degrad ation model has been 
hypothesized by the authors in earlier publications f. 5, 8 3 . The model allows a 
single parameter (denoted "characteristic degradation time") to describe the 
power degradation behavior of an electric propulsion spacecraft to a degree 
which fundamentally affects the solution to the trajectory optimization 
problem. In short, the power generated is degraded by a multiplicative 
damage factor q of the form 

-s/t 

q = e d, 

where r is the characteristic degradation time and s is the degradation time, 
d 

which is computed in the HILTOP program as the time integral of the density 
of damaging particles impinging on the solar arrays (for solar electric 
propulsion). The density of damaging particles is assumed to be a function of 
solar distance and array orientation. 

The characteristic degradation time is an engineering parameter that 
is determined experimentally. For example, by exposing a solar cell to the 
particle emission of a solar simulator and measuring the performance of the 
cell over a period of time, a reasonable value of can be estimated. 

The assumed exponential form of the degradation factor, although 
intended for use with SEP systems, is applicable for NEP systems as well. 

The principal difference is in the definition of s . The exponential form 
permits the evaluation of radio-isotope systems by defining s = 1 and letting 
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T, represent the time for the radioactivity to dissipate to l/e of its initial 
d 

level. A more complete exposition of this subject is given in »- 8 ] . 

Housekeeping Power . An option of simulating spacecraft house- 
keeping power has been added to the program. This option applies to solar 
electric propulsion with specified reference power. The housekeeping power 
is a specified constant power generated by the solar arrays and shunted away 
from the thruster power-conditioners and directly to the spacecraft payload 
for "housekeeping" purposes. The spacecraft model has been expanded by 
deleting the old, total propulsion system specific mass and replacing it with 
the specific mass of the solar arrays and the specific mass of the thruster 
subsystem. 

Declination Optimization . The program has been expanded to include 
the optimization of launch asymptote declination. The launch asymptote declina- 
tion optimization model was first hypothesized by the authors in the appendix of 
[l], and later a more thorough treatment of the subject was put forth in [9]. 

A solution to the problem of optimizing electric propulsion heliocentric 
trajectories, including the effects of geocentric launch asymptote declination on 
launch vehicle performance capability, hr - been developed using variational 
calculus techniques. The model of the launch vehicle performance includes a 
penalty associated with a non-easterly launch plus another penalty arising from a 
non-coplanar launch from the parking orbit. Provisions for range safety 
constraints are included. Optimal trajectories will generally have the launch 
excess velocity offset from the initial primer vector. The analysis describing 
the launch asymptote declination optim ization model is found in Section IV of 
this document. 

Swingby Continuation . Additional optional computations have been 
provided in which ballistic swingbys past the primary target may be simulated. 

In one mode of program operation, single swingbys past the primary 
target may be simulated to up to ten post-swlngby targets per case. This mode 
of operation was already in existence. 
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In another mode of program operation, multiple swingbys along a 
single trajectory may be simulated, first swinging past the primary target 
and then subsequently swinging past more targets downstream along the 
trajectory. One multiple swingby trajectory may be simulated per case. This 
mode of operation was recently added to the program. 

In either mode of operation, the following basic assumptions are made. 
The swingby continuation computations are independent of the trajectory leg 
leading up to the swingby tai'get, which may consist of an optimized electric 
propulsion trajectory segment (if the swingby target is the primary target), 
except that the arrival V and arrival time at the swingby target are used in the 
determination of the swingby passage conditions. Each swingby maneuver is 
calculated under the assumption of the patched-conic approximation, and the 
swingby planet's sphere-of-influence is assumed to have zero radius as seen 
from interplanetary space and infinite radius as seen from the planetary vantage 
point. The passage time in the swingby planet’s sphere-of-influence is neglected 
(taken to be zero in the heliocentric frame). 

Each swingby maneuver may be either unpowered or powered. The 
unpowered swingby solutions are embedded in the wider class of powered-swingby 
solutions, tending to appear in pairs which are separated by a region of braking 
powered swingbys. 

The powered swingby maneuver is restricted to occur at the mutual 
perifoci of the approach and departure hyperbolic arcs; the powered phase is 
impulsive and the thrust is collinear (pro or con) to the velocity at closest 
approach. Whether the swingby is powered or unpowered, the trajectory segment 
leading up to the swingby planet has been pre-determined, this being the method 
by which the program has been designed to obtain swingby solutions. Therefore 
the swingby time and the arrival hyperbolic excess velocity are known. 

A basic assumption of the powered swingby algorithm used by the 
program is that the flight time from the swingby planet to the next target is 
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specified. This being so, the program is able to converge, by iteration, on 
some ballistic trajectory from the swingby planet to the next target having 
the specified transfer time, implying that the departure hyperbolic excess 
velocity at the swingby planet is thereby determined. Therefore, the helio- 
centric trajectory before and after the swingby planet is determined, and it 
then remains to perform the required computations pertaining to the hyperbolic 
arcs within the swingby planet's sphere of influence. The computations are 
outlined in a companion document [lO] which is being published concurrently 
with this report. 

The unpowered swingby maneuver is considered to be a powered 

swingby having AV= 0. The program adjusts the post-swingby heliocentric 

trajectory segment, by iteration, until the swingby departure magnitude 

equals the given arrival V magnitude. The primary independent variable in 

00 

this iteration is the post-swingby transfer time to the specified target, which 
was held constant in the powered swingby case. 

The program can generate multiple-revolution ballistic arcs, and a 
particular solution obtained by the program may not be unique, even for the 
same transfer time. All solutions are reachable, however, by means of in- 
putting an appropriate initial velocity guess for the trajectory segment in 
question. 

Deep Space Burn . In simulations of trajectories which are all -ballistic, 
the program is now capable of simulating a single deep space burn, or impulsive 
velocity -change, at any point prior to arrival at the primary target. The three 
components of the incremental velocity AV are independent variables of the 
boundary value problem, such that, at a specified time, the spacecraft velocity 
is Incremented: 

R + = R‘ + 

The optional existence of a deep space burn provides greater targeting 
flexibility in simulations of multiple target missions. 
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V Optimization in LVI Mode . The optimization of the launch excess 

"?0 

velocity V o when using the Launch Vehicle Independent (LVI) mode of simulation 
Is accomplished when the initial primer vector is forced to vanish: 

A = 0. 

o 

This is accomplished by setting the values of these three independent variables of 
the boundary value problem to zero and turning their triggers off; the three com- 
ponents of the departure heliocentric velocity become independent variables in- 
stead. The program has been augmented with special logic to circumvent the 
numerical singularity associated with the null primer vector. 

Print Expansion . The standard trajectory block print has been expanded 
to include target-relative position and velocity coordinates of the spacecraft, 
and comet nuclear and total magnitudes as seen by both the spacecraft and 
Earth. A more detailed description is given in Cio] . 

Extra-Ecliptic Missions . Extra-ecliptic mission simulations now 
involve launches from the Earth in which the Earth’s ephemeris is generated 
by the program's analytic ephemeris capability; previously, extra-ecliptic 
missions were generated simply by starting the trajectory on the x-axis at one 
astronomical unit from the sun. This improved extra-ecliptic capability 
allows the launch date to be optimized together with the launch asymptote 
declination. Also, an additional set of boundary conditions has been added 
to the program for simulating extra-ecliptic missions, and these are given 
in [10: . 

Miscellaneous Improvements . Ephemerides for several additional comets 
and asteroids have been included in the program's analytic ephemeris capability, 
increasing the total number of possible targets to fifty-one. The ephemeris 
arrays have been expanded to allow up to seventy targets. 

The program option of generating a ballistic trajectory as an initial 
guess for an electric propulsion mission has been expanded to allow multiple- 
revolution trajectories. 
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The capability has been added of monitoring the proximity of any 
given spacecraft trajectory to any object in the solar system. 

Another capability has been added of specifying enforced coast 
phases for selected intervals throughout a given mission. 

Finally, a major size-reduction of the HILTOP program was 
completed, from 436K to 326K hexadecimal locations, a reduction of 25%, 
allowing faster turn-around time on the GSFC IBM 360/91. 
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VI. 


QUICKTOP 1II/CHEBYTOP III COMPUTER PROGRAM USAGE 


[6] [7] 

The QUICKTOP III and CHEBYTOP III computer programs have 
been obtained from the NASA Ames Research Center, Moffett Field, California, 
and converted to run on the IBM 3G0/91 computer at the Goddard Space Flight 
Center. These trajectory optimization computer programs were used initially 
In the Sert III comet rendezvous mission study discussed in Section II, and 
aided in establishing the feasibility of using the Sert III spacecraft for these 
comet rendezvous missions. However, it was determined early in the study 
that the QUICKTOP/CHEBYTOP programs were consuming as much or more 
machine time as the HILTOP^ program. For this reason, a decision was 
made to solely use the HILTOP program for the comet rendezvous study, since 
the authors are more familiar with this computer program. The QUICKTOP/ 
CHEBYTOP computer programs remain ready for use at the Goddard Space 
Flight Center by any interested persons. 
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VII . ERROR ANALYSIS OF HIGH-THRUST MANEUVFRS 

High thrust burns can result in large postburn errors when applied to 
spin-stabilized spacecraft. In addition to the simulation of the burn, an 
algorithm predicts the expected errors after high thrust maneuvers. The pro- 
cedure incorporates navigational uncertainties in the pre-burn state, and attitude 
and magnitude errors during the burn to estimate errors in the resulting orbital 
elements. The work was performed by W. Bjorkman. 

To add to the concreteness of the analytical development, the Synchronous 
Meteorological Satellite (SMS) mission is discussed as a sample case. The apogee 
burn of the SMS mission is performed by firing the solid rocket of a spin-stabilized 
spacecraft near apogee of the transfer orbit. The rocket burns for about 23 
seconds, imparting a delta-velocity of about 1723 m/sec along the direction of 
the spin axis. The time of firing the rocket and the attitude of the spin axis are 
controllable. Appropriate values for these control parameters (firing time and 
attitude angles) can be determined from a scanning procedure. Any procedure 
used to establish control parameters will use "best estimates" of the transfer 
orbit state and the expected delta-velocity from the rocket. The estimated trans- 
fer orbit state will always differ from the true orbit state because of navigation 
errors (i.e., errors caused by neglected error sources, incomplete modeling or 
measurement errors in the orbit determination process). The solid rocket will 
not deliver exactly the delta-velocity stated by the manufacturer. In addition to 
these errors, the knowledge of spin-axis attitude will be in error by uncertainties 
in the attitude determination process. An algorithm was developed for assessing 
the effects of navigation errors and burn errors (i.e. , delta-velocity and attitude 
determination errors) on the achievement of mission objectives. Attitude control 
errors are not considered, nor are errors in firing time. The algorithm makes 
use of an impulsive burn model and linear propagation of errors through the 
burn. A particular set of mission objectives is assumed without assuming a 
linear relationship between mission constraint parameter variations and naviga- 
tion/burn errors. The validity of the linear error propagation assumption is 
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demonstrated with a numerical example. 

Navigation errors are specified by means of a 6x6 covariance matrix, 

P , referred to th anchor epoch, t , and a suitable coordinate frame. P 
o o o 

should include uncertainties due to all orbit determination errors, and not simply 
those uncertainties caused by measurement noise. The Cartesian state error, 
x q , is the difference between the true state at epoch and the "anchor vector". 

We assume 


and 


e<x ) = o 

o 


£(x x T ) = P 
o o o 


( 1 ) 

(2) 


where £ is the expectation operator. References to equations in this section 
of the report apply to equations in this section only. Small deviations in the 
state propagate linearly along the transfer orbit between t and the firing time, 
tj, by means of the "mean conic" state ti'ansition matrix, <P(tj.; t Q ). 


l o )x ~o 


(3) 


This transition matrix is evaluated from the anchor vector and the osculating 
state (integrated from the anchor vector) at t^., using the average reciprocal 
semi-major axis at the two terminals in calculation of the incremental regulariz- 
ing anomaly. Table VII -1 illustrates numerically the adequacy of equation 
(3) for propagation of errors along the SMS-A transfer orbit. Linearity errors 
may be smaller than those shown in the table, because a coordinate system in- 
consistency was detected after the tabulated data were generated. The transition 
matrix is deterministic and therefore commutes with the expectation operator. 


e(px Q ) =v5£(x q ) = o 

(4) 

P f = e (x f x f T ) =P e <V> T =<f>P o <p T 

(5) 


The effect of navigation errors on the state at firing time Is thus given by P^ in 
equation (5). 
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Apogee burn errors are specified as "proportional” and "pointing" errors. 
The proportional error is a delta-velocity error which is proportional to the ex- 
pected total delta-velocity, AV. Random behavior is assigned to the proportional 
error. 


£(k ) = 0 

V 

(6) 

2 2 . — .9 

e(k ) =• a AV 

V V 

(7) 


Attitude errors contribute a delta-velocity error approximately in the plane 
normal to AV. Unit basis vectors in that plane are defined by 


E = 


and 


N = 


2 

_ AV 


6v i 
V 0/ 


x E 


|a v I 


( 8 ) 


0) 


Pointing error is a circular error probability number assigned to attitude de- 
termination accuracy. It is treated as having a normal distribution about zero 
mean in each of two orthogonal directions. 


e(k ) = e (k ) = 0 
e n 

e(k 2 ) = £(k 2 ) = a 2 |Av | 2 
e n a 

The delta -velocity vector (with errors included) may be written 

r AV+k E+k N - 

AV = f|Avj ) r _ « .. . 1 

W |AV + k E+k N| J 
e n 1 


(10) 

(ID 


(12) 
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To linearize the burn we write the delta-velocity error In the form 

AV= M I r | = Mr 


V 


n 


( 13 ) 


where r , r and r are uncorrelated random numbers of zero mean and unit 
v e n 

variance, defining 




|av| 


/o rN 
v v 

0 r 
a e 





Vo r , 
'an 


( 14 ) 


The covariance matrix of delta-velocity errors Is 

** rn rp rp rp 

C(AV AV ) = M£(fr )M =MM (15) 

where M is defined by the following equation 


M = 


5 AV 


d(r ,r ,r ) 
' v e n 


r =r =r =0 
v e n 



a6) 


The adequacy of the linear model for burn errors is illustrated numerically in 
Table VII -2. The comparisons were obtained by multiplying random samples 
of r by M of equation (16) and differencing these with (AV of equation (12) 

- AV). The input standard deviations were: O y = .0025 and = (.43/57.3)°. 
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Table VII -2 

Linear Propagation Verification: Durn Errors 

XA - M r __ 

XE = XA - (AV - AV) 

(In meters/second) 


Sample 

e 

X 

e 

y 

e 

z 

1 XA 

. 3336 

-.0109 

-.3015 

XE 

-.0000 

-.0001 

.0000 

2 XA 

-.8110 

1.5135 

.7079 

XE 

-.0001 

- .0001 

-.0000 

3 XA 

2,6601 

-7.2964 

-3.6459 

XE 

.0022 

.0004 

.0010 


The post-burn state Is computed impulsively by adding AV to the pre- 
burn velocity. 


X =X f + 



(17) 


Errors are added the same way, 



(18) 


The covariance matrix of post-burn errors, P + , is then computed adding 

** ^ rp 

C(AV AV ) of equation (15) to of equation (5). Navigation and burn errors 
are assumed uncorrelated. 


+ 

P = P f + 


0 ° T 1 

0 MM 1 _ 


(19) 


A vector, of mission constraint parameters was defined for the SMS 
mission. Each element of 0 can be computed unambiguously from X + . The 
parameters are: 
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required three-impulse trim velocity required to circularize 
at synchronous radius with a specified inclination and node 


lb 

* 3 

*4 


longitudinal drift rate 
eccentricity of the pre-trim orbit 
inclination of the pre-trim orbit 
ascending node of the pre-trim orbit 


Table VII -3 compares constraint vector errors computed from linearly- 
propagated errors with errors propagated by "non-linear" means. 


Table VII -3 


Linear Propagation Verification: Constraint Parameters 
PSI = ^(X + -* <px' Q + ( M - ) ) linear 




ERR = XA 

- (X + ) 

linear 

$ error 


Sample 


TRIMV 

DRIFT 

e 

i 

O 

1 

PSI 

253.59 

64.97 

.057456 

.659 

151.932 


ERR 

.06 

.02 

.000013 

.000 

.008 

2 

PSI 

284.83 

66.78 

.045527 

1.117 

136.212 


ERR 

.06 

.02 

.000020 

.000 

.003 

3 

PSI 

317,76 

68.99 

.041745 

1.G1G 

13C.577 


ERR 

.07 

.02 

.000026 

.000 

.001 

4 

PSI 

252.05 

64.81 

.058739 

.645 

153.628 


ERR 

.13 

.04 

.000027 

.000 

- .020 

5 

PSI 

283.13 

66.59 

.046765 

1.097 

136.928 


ERR 

.13 

.04 

.000041 

.000 

.006 

6 

PSI 

315.89 

68.77 

. 042661 

1.592 

130.961 


ERR 

.14 

.04 

.000052 

.000 

- .002 
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Variation of some of these parameters with x is very non-linear, so a 
Monte Carlo procedure was implemented as part of the error analysis algorithm. 
The valid assumption of linearity of error propagation from anchor epoch through 
the burn (equations (3), (13), and (18)) makes is reasonable to sample post-burn 
errors directly from P of equation (19). 

Statistical characteristics of the constraint parameter errors are com- 
puted from the Monte Carlo samples. These characteristics (as implemented 
to test the algorithm) a: 


1. 

minimum value 


2. 

maximum value 


3. 

mean 


4. 

standard dev’ation 


5. 

probability that | $ | 

s£(£ input range) 

6. 

probability that 1 0 | 

s 2C 

7. 

probability that 1 0 | 

£ 36 


A frequency histogram of 20 equal algorithm intervals about zero was 
also implemented for the test. All tests indicate that the algorithm as described 
is a valid one for error analysis of high thrust maneuvers. 

The algorithm for error analysis of high-thrust maneuvers may thus be 

stated: 

1. Propagate the navigation error covariance matrix to firing time 
(equation (5)). 

2. Add the velocity-impulse error covariance matrix, M, to the 
navigation error covariance matrix at firing time (equation (19)). 

3. Sample post-ourn state errors from the resulting state error 
covariance matrix using a random number generator. 

4. Add the sampled state errors to the errorless post-burn state. 
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5. Compute the mission constraint parameter error by differencing 
the value computed with the erroneous state ai.d that computed 
with the errorless state. 

6. Accumulate statistics of constraint parameter errors, 

7. Return tc step 3 a specified number of times. 

8. Display statistics of constraint parameter errors. 
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APPENDIX A 


This appendix consists of detailed tabular data pertaining to Section II 
on comet rendezvous missions. A glossary of the tabulated parameters precedes 
the tables, which have headings sufficient to define their contents. 

The explanations and comments below may be helpful in clarifying the 
definitions and usefulness of some of the tabulated quantities. Refer also to the 
glossary mentioned in the preceding paragraph. These comments are in no parti- 
cular order. 

(1) The quantities r and 0 may be used to determine (on all but 

' com com 

the first table "Launch Date Variation") the rendezvous communication 
distance and angle along the 
orbit of the comet, i.e. , after 
the rendezvous is accomplished 
and while the spacecraft is travel- 
ing along with the comet. Also, 
the pj. tabulation effectively dis- 
plays the power available to the 
spacecraft as it moves along with the comet. 

(2) The first table, labelled "Launch Date Variation", may be used to obtain 
a rough idea of the mission window width for the Giacobini-Zinner/1985 
rendezvous mission using a Titan III E/Centaur launch vehicle; the re- 
sults of a mission window analysis would depend on the analyst’s criteria. 

(3) In the burn time, t^, tabulation, "cont. " means "continuous" which stands 
for continuous thrusting throughout the pre-rendezvous portion of the 
mission; t^ = tj. See the discussion in Section II of the main report con- 
cerned with single-engine throttling for the proper method of computing 
individual engine on-time requirements. 


SUN 



SPACECRAFT 
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(4) Pertaining to 0 . and 0 , the thrust cone angle 0 is depicted by 

' ' min max 


the following sketch: THRUST 



m = m + m + m + m 
o ps p t net 


where m pg is the propulsion system mass (see main text, Section II, 

for discussion), equal to 262.6 kilograms, m p and m^ are tabulated, 

and m = 0.1 x m . m was computed by the HILTOP program as a 
t p o 

function of the launch characteristic speed v c according to the expression: 


m 

o 


, ^“ V /bJ 

= b^e c 2 



For m in kilograms and v in meters per second, the coefficients 
o c 

are: 

b b b 

12 3 

Titan III E/Centaur 167238.95 3480.2038 1753.6965 

Shuttle/Transtage 2859382.94 1715.7632 1199.9231 

(6) The importance of the launch asymptote declination 6 and parking orbit 

inclination i are discussed in Section IV of the main text. When 

|6| £ 28°. 5, the launch is due-east, and, for possible comparison with 

other studies, the tabulated m and v columns give the dependence of 

o 00 

launch vehicle injected mass m Q on departure excess speed v o> This 
is not true when |6 | > 28°. 5 because then other factors contribute to 
the launch phase characteristic speed. 
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GLOSSARY OF TABULATED PARAMETERS 


AD 

r mln 

r 

max 

r 

com 

6 

com 

0 . 

mm 

0 

i ax 

6 

i 

LD 


m 


m 


net 


Arrival date, days from perihelion passage. 

Minimum solar distance encountered along trajectory, AU. 
Maximum solar distance encountered along trajectory, AU. 
Earth-spacecraP distance at arrival, AU. 

Sun-Earth-spacecraft configuration angle at arrival, degrees. 
Minimum thrust cone angle required, degrees. 

Maximum thrust cone angle required, degrees. 

Change in longitude during the mission, degrees. 

Launch asymptote declination, degrees. 

Inclination of launch parking orbit, degrees. 

Launch date, calendar date. 

Mission duration, days. 

Total propulsion system on time, days. 

Initial spacecraft mass, kilograms. 

Low thrust propellant mass, kilograms. 

Net spacecraft mass, kilograms. 

Power developed by arrays at target arrival, kilowatts. 
Launch hyperbolic excess speed, meters per second. 


v 




G(acobtnI-Z(nner(85) Rendezvous Titan in E/Centaur/Sert in Launch Date Variation 


AD 

r m!n 

r 

max 

r 

com 

e 

com 

0 f 
min 

0 

max 

AX 

6 

l 

LD 


*b 

m 

o 

m 

P 

m 

net 

P t 

tr 

m 

0 

• 844 

3.992 

.45 

79.6 

50.9 

159.3 

470.7 

44.27 

32.50 

6/8.5/81 

1550 

1394.5 

1127.7 

469.5 

348.6 

8.37 

8053.7 

0 

.834 

4.045 

.45 

79.6 

49.5 

163.7 

461.2 

38.44 

32.50 

6/18.5/81 

1540 

1406.9 

1134.0 

464.9 

360.0 

6.37 

8619.6 

0 

.840 

4.073 

.45 

79.6 

46.1 

166. 7 

451.6 

31.67 

30.98 

6/28.5/81 

1530 

1417.7 

1135.5 

459.9 

367.0 

8.37 

8813.5 

0 

.855 

4.085 

.45 

79.6 

41.1 

16G.1 

442.1 

23.93 

28.50 

7/8.5/81 

1520 

1424.1 

1137.5 

453.5 

375.9 

8.37 

8829.7 

0 

.878 

4.082 

.45 

79.6 

35.2 

160.2 

432.5 

15.10 

23.50 

7/J8.5/81 

1510 

1430.2 

1147.8 

448.7 

391.6 

8.37 

8809. 8 

cr> 0 

.905 

4.0G7 

.45 

79.6 

29.2 

151.0 

423.0 

5,24 

28.50 

7/28.5/81 

1500 

1433.7 

1167.5 

444.1 

416.4 

8.37 

8772.0 

to o 

.935 

4.041 

.45 

79.6 

24.0 

139.6 

413.4 

-5.32 

28.50 

8/7.5/81 

1490 

1434.5 

1198.8 

439.7 

452.6 

8.37 

8712.4 

0 

. 904 

4.011 

.45 

79.6 

20.9 

127.4 

403. 8 

-15.26 

29.50 

8/17.5/81 

1480 

1432.6 

1242.5 

435.1 

501.4 

8.37 

8629.6 

0 

.9K5 

3.990 

.45 

79.6 

20.3 

116.0 

394.2 

-22.48 

28.50 

8/27.5/81 

1470 

1427.5 

1295.1 

429.5 

560.0 

6.37 

6331.0 

0 

. 996 

3.9S4 

.45 

79.6 

21.0 

106.7 

334.5 

-26.17 

28.50 

9/6.5/81 

1460 

1418.9 

1317.0 

422.5 

619.6 

8.37 

8* 34. 7 

0 

1.001 

3.992 

.45 

79.6 

21.7 

101.9 

374.8 

-27.16 

28.50 

9/16.5/81 

1450 

1407.3 

1385.8 

414.2 

667.6 

8.37 

8363.3 

0 

1.001 

4.009 

.45 

70.6 

21.8 

103.5 

365,0 

-26.42 

28.50 

9/26.5/81 

1140 

1393.6 

1100.5 

405.5 

691.8 

8.37 

8336.5 

0 

1.000 

4.026 

.45 

79.6 

21.4 

101.4 

355.2 

-24.50 

28.50 

10/6.5/81 

1430 

1379.6 

1584.8 

397.8 

684.6 

8.37 

S365.3 

0 

.997 

4.045 

.45 

79.6 

20.4 

99.3 

345.3 

-21.52 

28.50 

10/16.5/81 

1420 

1367.6 

1339.1 

391.7 

645.6 

8.37 

8149.4 

0 

.994 

4.065 

.45 

79.6 

19.1 

97.5 

335.3 

-17.32 

28.50 

10/26.5/81 

1410 

1358.1 

1269.3 

386. 4 

581.6 

8.37 

8579.3 

0 

.991 

4.098 

.45 

79.6 

17.7 

96.1 

325.3 

-11.46 

29.50 

11/5.5/81 

1400 

1351.2 

1185.5 

381.1 

503.7 

8.37 

8737.7 

0 

.989 

4.114 

.45 

79.6 

16.8 

95.4 

315.3 

-3.75 

28.50 

11/15.5/81 

1390 

1348.4 

1098.8 

376.3 

422.2 

8.37 

8904 5 

0 

.987 

4.149 

.45 

79.6 

16.8 

96.0 

305.2 

5.04 

28.50 

11/25.5/81 

1380 

1353.9 

1017.1 

373.9 

343.2 

8.37 

9064.5 
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4-Year Missions 


AO 

r mln 

r 

max 

r 

com 

Q 

com 

^mln 

0 

max 

AX 

6 

I 

LD 


*b 

m 

o 

m 

P 

m 

net 

P t 

V 

m 

0 

1.001 

4.013 

.45 

79.6 

21.7 

102.9 

362.2 

-26.00 

28.50 

9/29.3/81 

1437. 2 

1389.6 

1399.2 

403.1 

693.1 

8.37 

8338.8 

-10 

1.001 

4.013 

.47 

80.0 

23.0 

105.8 

351.6 

-26.54 

28.50 

9/29.3/81 

1427.2 

1385.1 

1384.3 

396.2 

685.9 

8.26 

836G.1 

-20 

1.001 

4.017 

.51 

81.8 

24.7 

109.3 

341.1 

-27.07 

28.50 

9/29.2/81 

1417.3 

1383.2 

1371.5 

390.9 

679.0 

7.97 

8389.6 

-30 

1.001 

4.020 

.57 

84.3 

26.7 

113.4 

331.0 

-27.56 

28.50 

9/28.9/81 

1407.6 

1384.0 

1361.2 

387.4 

672.4 

7.53 

8408.6 

-40 

1.000 

4.022 

.64 

87.2 

28.9 

118.0 

321.3 

-27.98 

28.50 

9/28.5/81 

1398.0 

1387.9 

1353.4 

385.7 

666.5 

6.99 

8423.0 

-50 

1.000 

4.026 

.71 

90.1 

30.5 

122.4 

312.2 

-28.35 

28.50 

9/28.2/81 

1388.3 

cont. 

1344.1 

382.1 

661.1 

6.43 

8440.2 

-60 

1.000 

4.042 

.79 

92.8 

29.8 

126.9 

302.7 

-28.77 

28.50 

9/28.8/81 

1377.7 

cont. 

1324.4 

372.8 

651.8 

5.86 

8476.1 

-70 

.998 

4.067 

.88 

95.2 

27,2 

130.7 

292.7 

-29.13 

28.50 

9/30.8/81 

1365.7 

cont. 

1294.7 

362.9 

632.9 

5.33 

8529.7 

-80 

.997 

4.099 

.98 

96.9 

23.5 

134.1 

283.1 

-29.57 

28.93 

10/3.2/81 

1353.3 

cont. 

1254.3 

353.4 

602.9 

4.85 

8G02.2 

-90 

.995 

4.133 

1.08 

98.0 

19.5 

137.1 

274.3 

-30.10 

29.44 

10/5.5/81 

1341.0 

cont. 

1205.8 

344.7 

564.1 

4.41 

8690.1 

-too 

.994 

4.167 

1.19 

98.2 

15.7 

139.6 

266.6 

-30.64 

29.97 

10/7.4/81 

1329.1 

cont. 

1153.0 

336.8 

520.0 

4.03 

8787.1 

■110 

.991 

4.201 

1.32 

97.6 

12.1 

141.8 

259.9 

-31.12 

30.43 

10/8. 9/81 

1317.6 

cont. 

1099.2 

329.7 

474.0 

3.G8 

8887.7 

■120 

.989 

4.234 

1.45 

96.0 

9.0 

143.6 

254.2 

-31.48 

30.78 

10/10.0/81 

1306.5 

cont. 

1047.1 

323.5 

428.7 

3.38 

8937.1 

■130 

.985 

4.265 

1.60 

93.7 

6.1 

145.2 

249.4 

-31.69 

30.98 

10/10.7/81 

1295.9 

cont. 

996.8 

318.0 

384.4 

3.12 

9085.1 

*140 

.982 

4.294 

1.76 

90.7 

3.4 

146.6 

245,3 

-31.71 

30.99 

10/11.0/81 

1285.5 

cont. 

950.2 

313.1 

343.1 

2.88 

9178.2 

■150 

.978 

4.321 

1.93 

87.0 

0.9 

147.8 

241.9 

-31.53 

30.81 

10/11.0/81 

1275.5 

cont. 

907.0 

309.0 

304.6 

2.67 

9266.7 


3-Year Missions 


0 

.999 

3.135 

.48 

79.3 

30.7 

123.4 

360.9 

-34*88 

32.50 

-1C 

.998 

3.143 

.48 

79.8 

32.6 

127.0 

350.8 

-35.29 

32.50 

-20 

.996 

3.159 

.52 

81.5 

34.2 

130.4 

340.6 

-35.65 

32.50 

-30 

.992 

3.190 

.57 

84.0 

34.0 

133.6 

329.9 

-35.88 

32. 50 

-40 

.985 

3.235 

.64 

86.9 

31.1 

136.8 

318.7 

-35.87 

32. 50 

-50 

.977 

3.285 

.72 

89.8 

26.7 

139.7 

308.2 

-35.61 

32.50 

-60 

.967 

3.333 

.80 

92.5 

22.2 

142.2 

298.9 

-35,11 

32. 50 

-70 

.955 

3.378 

.89 

94.8 

18.1 

144.4 

291.1 

-34.41 

32.50 

-80 

.943 

3.420 

.98 

96.6 

14.3 

146. 2 

284.5 

-33.55 

32.50 


10/1.4/82 

1070.1 

cont. 

1336.7 

434.4 

596.3 

8.37 

8392.9 

9/30. 9/82 

1060.6 

cont. 

1306.9 

423.0 

579.0 

8.26 

b436.9 

9/30. 5/82 

1051.0 

cont. 

1272.4 

411.7 

557.0 

7.97 

8489.9 

9/30.8/82 

1040.7 

cont. 

1227.5 

399.9 

525.0 

7.52 

8567.2 

10/1.9/82 

1029.6 

cont. 

1167.9 

387.9 

478.6 

6.99 

8682.4 

10/2.9/82 

1018.6 

cont. 

1098.8 

376.6 

422.0 

6.43 

8826.1 

10/3.4/82 

1008.1 

cont. 

1028.0 

566.4 

362.3 

5.86 

8980.7 

10/3.2/82 

998.3 

cont. 

960.6 

357.5 

304.8 

5.33 

9131.5 

10/2.4/82 

989.0 

cont. 

80S. 6 

349.7 

251.4 

4.85 

9260.6 
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Titan III E/Centaur/Sert m 



Borrelly(87) Rendezvous 


AD 

r rato 

r 

max 

r 

com 

6 

com 

*mln 

0 

max 


6 

1 

ID 

0 

.978 

4.026 

.51 

127.5 

36.7 

90.0 

360.0 

37.83 

32.50 

12/1.8/83 

-10 

.978 

4.029 

.49 

131.8 

36.6 

90.1 

353.3 

37.87 

32.50 

1 2/1.6/83 

-20 

.978 

4.031 

.50 

133.2 

36.7 

90.1 

346.6 

37.92 

32.50 

12/1.3/83 

-30 

.977 

4.034 

.53 

132.1 

36.9 

90.3 

340.1 

37.96 

32.50 

11/30. 9/83 

-40 

.977 

4.036 

.58 

129.5 

37.2 

90.9 

333.8 

38.01 

32.50 

11/30.4/83 

-50 

.976 

4.039 

.65 

126.6 

37.5 

93.7 

327.6 

38.05 

32.50 

11/29.9/83 

-60 

.970 

4.042 

.72 

124.0 

37.9 

97.7 

321.6 

38.08 

32.50 

11/29.3/83 

-70 

.975 

4.046 

.80 

121.6 

38.0 

102.6 

315.8 

38.11 

32.50 

11/28.9/83 

-80 

.975 

4.052 

.89 

119.4 

37.6 

108.3 

309.9 

38.12 

32.50 

11/28.7/83 

-90 

.974 

4.062 

.98 

117.3 

36.0 

114.4 

303.7 

38.08 

32.50 

11/29.1/83 

-100 

.974 

4.079 

1.08 

115.1 

32.8 

119.8 

296.8 

37.95 

32.50 

11/30.5/83 

-110 

.974 

4.101 

1.18 

112.7 

28.5 

124.3 

289.8 

37.71 

32.50 

12/2.3/83 

-120 

.974 

4.126 

1.29 

110.1 

24.0 

128.0 

282.9 

37.39 

32.50 

12/-S.3/83 

-130 

.973 

4.152 

1.41 

107.2 

19.9 

131.1 

276.6 

37.00 

32.50 

12/6.1/83 

-140 

.972 

4.178 

1.54 

104.0 

16.4 

133.6 

270.8 

36.55 

32.50 

12/7.6/83 

-150 

.971 

4.204 

1.67 

100.4 

13.5 

135.9 

265.6 

36.05 

32.50 

12/8.8/83 


0 

.969 

3.180 

.51 

128.0 

40.5 

102.3 

360.8 

38.87 

32.50 

11/30.2/84 

-10 

.968 

3.186 

.49 

132.3 

41.5 

105.8 

354.4 

38.84 

32.50 

11/29. 7/84 

-20 

.966 

3.193 

.50 

133.6 

42.5 

110.0 

348.0 

38.79 

32.50 

11/29.2/84 

-30 

.965 

3.203 

.53 

132.3 

43.3 

114.9 

341.6 

38.68 

32.50 

11/28.7/84 

-40 

.963 

3.221 

.58 

129.6 

43.6 

120.4 

334.8 

38.43 

32.50 

11/29.6/84 

-50 

.960 

3.249 

.65 

126.6 

42.3 

125.9 

327.1 

37.87 

32.50 

11/29.6/B4 

-60 

.957 

3.265 

.72 

123.9 

39.9 

130.5 

318.9 

37.00 

32.50 

12/1.2/84 

-70 

.953 

3.321 

.80 

121.5 

37.0 

134.3 

311.2 

35.90 

32.50 

12/2.6/84 

-80 

.948 

3.356 

.89 

119.3 

34.2 

137.4 

304.3 

34.68 

32.50 

12/3.4/84 


4- Year Missions 



S 

m 

o 

m 

P 

m 

net 

P t 

▼ 

m 

1477.7 

cont. 

1386.7 

423.7 

658.1 

5, 62 

8167.1 

1467.9 

cont. 

1371.6 

416.2 

651.2 

5.59 

8192.9 

1458.2 

cont. 

1356.9 

408.8 

644.6 

5.50 

6218.0 

1448.6 

cont. 

1342.9 

401.7 

638.4 

5.35 

8242.0 

1439.1 

cont. 

1329.6 

394.9 

632.6 

5.16 

8264.9 

1429.6 

cont. 

1317.0 

388.3 

627.2 

4.94 

8286.6 

1420.2 

cont. 

1304.8 

382.1 

621.9 

4.69 

8307.9 

1410.6 

cont. 

1292.4 

376.0 

616. 2 

4.44 

8329.9 

1400.8 

cont. 

1278.8 

369.9 

609.3 

4.19 

8355.6 

1390.4 

cont. 

1261.6 

363.4 

599.3 

3. 94 

8391.5 

1379.0 

cont. 

1238.1 

356.4 

583.5 

3.70 

8445.7 

1367.2 

cent. 

1207.2 

349.1 

560.6 

3.47 

8520.1 

1355.2 

cont. 

1170.2 

342.1 

531.3 

3.26 

8610.3 

1343.4 

cont. 

1129.5 

335.6 

497.7 

3.06 

8710.3 

1331.9 

cont. 

1087.0 

329,6 

461.8 

2.88 

8814.6 

1320.7 

cont. 

1044.5 

324.1 

425.3 

2.71 

8919.7 


3-Year Missions 


1113.1 

cont. 

1286.0 

425.6 

555.3 

5.62 

8290.3 

1103.6 

cont. 

1265.0 

418.0 

542.6 

5.59 

8333.4 

1094.1 

cont. 

1242.5 

410.5 

528.4 

5.50 

8381.4 

1084.6 

cont. 

1217.3 

402.9 

511.5 

5.35 

8439.7 

1074.6 

cont. 

1186.6 

394.9 

489.6 

5.16 

8516.7 

1063.6 

cont. 

1146.5 

395.9 

459.5 

4.94 

8630.1 

1052.0 

cont. 

1097.3 

376.5 

420.6 

4.69 

8773.7 

1010.6 

cont . 

1043.7 

367.7 

376.6 

4.44 

8925.3 

1029.8 

cont. 

989.7 

359.8 

331.3 

4.19 

9067.8 
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Tempel 11(88) Rendezvous Titan It! E/Centaur/Sert ID 


AD 

r mfn 

r 

max 

r 

com 

0 

com 

0 

mtn 

0 

max 

AX 

6 

I 

LD 

0 

1.010 

3.919 

.96 

89.4 

64.1 

98.0 

371.0 

-5.39 

28.50 

7/21.5/84 

•10 

1.015 

3.923 

.92 

91.9 

64.1 

99.2 

363.1 

-5.33 

28.50 

7/22.1/84 

•20 

1.015 

3.927 

.89 

94.9 

64.1 

100.7 

355.3 

-5.68 

28.50 

7/22.7/84 

•30 

1.015 

3. 931 

.86 

98.6 

84.1 

102.8 

347,6 

-5.82 

28.50 

7/23.3/84 

-40 

1.015 

3.936 

.83 

103.0 

64.1 

104.8 

340.2 

-5.94 

28.50 

7/23.8/84 

•50 

1.015 

3.941 

.81 

108.2 

64.1 

107.3 

333.1 

-6.04 

28.50 

7/24.3/84 

-60 

1.015 

3.916 

.79 

114.3 

64.1 

110.0 

326.4 

-6.10 

28.50 

7/24.9/84 

•70 

•80 

•90 

1.015 

3. 952 

.78 

121.2 

64.2 

112.8 

319.9 

-6.11 

28.50 

7/25.4/84 

1.015 

3. 963 

.78 

137.3 

64.5 

119.5 

308.3 

-6.37 

28.50 

7/26.0/84 

-100 

1.015 

3.968 

.79 

145.5 

64.8 

123.1 

303.2 

-6.73 

28.50 

7/26.1/81 

-no 

1.014 

3.976 

.83 

152.4 

65.1 

126.4 

208.3 

-7.07 

28.50 

7/26.3/84 

•120 

1.014 

3.987 

.88 

156.1 

65.0 

129.4 

293.6 

-7.43 

28.50 

7/26. G/84 

•130 

1.014 

4.003 

.95 

155.0 

64.2 

131.6 

288.8 

-7.87 

28.50 

7/27.3/84 

-uo 

1.014 

4.025 

1.04 

149.8 

61.7 

133.3 

283.9 

-8.51 

28.50 

7/28.4/84 

-150 

1.013 

4.045 

1.14 

142.7 

57.7 

135.0 

279.0 

-9.26 

28.50 

7/29.6/84 

-160 

1.013 

4.063 

1.27 

134.9 

52.8 

136.7 

274.5 

-9.91 

28.50 

7/30.8/84 

-170 

1.012 

4.079 

1.41 

127.1 

48.1 

138.4 

270.6 

-10.28 

28.50 

7/31.5/84 

-180 

1.010 

4.093 

1.56 

119.4 

43.6 

140.0 

267.2 

-10.38 

28.50 

7/31.8/84 

•100 

1.009 

4.106 

1.73 

111.9 

39.2 

141.7 

264.3 

-10.24 

28.50 

7/31 . 9/84 

-200 

1.007 

4.118 

1.90 

104.7 

35.0 

143.4 

261.8 

-9.91 

28.50 

7/31.6/84 


4- Year Missions 


t # 

t 

m 

m 

m 

p 

▼ 

f 

D 

o 

P 

net 

t 

m 

1518.0 

1387.1 

2361.7 

395.9 

1663.6 

5.45 

6690.9 

1507.4 

1387.3 

2353.1 

391.9 

1659.5 

5.43 

6704.9 

1496.8 

1388.2 

2344.1 

388.1 

1654.6 

5.35 

6719.6 

1486.2 

1390.1 

2334.8 

394.7 

1649.1 

5.23 

6734.8 

1475.7 

1393.2 

2325.3 

381.6 

1642.9 

5.08 

6750.3 

14G5.2 

1398.1 

2315.8 

379.0 

1636.3 

4.89 

6765.9 

1454.6 

1405.2 

2306.3 

376.8 

1629.2 

4.69 

6781.3 

1444.1 

1415.2 

2297.1 

375.3 

1621.7 

4.47 

6796.4 

1423.5 

cont. 

2274.2 

369.4 

1605.3 

4.03 

6834.0 

1413.4 

cont. 

2257.5 

363.7 

1594.8 

3.81 

6861.4 

1403.2 

cont. 

2238.1 

358.3 

1581.4 

3.60 

6893.2 

1392.9 

cont. 

2213.7 

352.9 

1562.9 

3.40 

6933.4 

1382.2 

cont. 

2180.0 

347.4 

1535.3 

3.22 

6988.9 

1371.1 

cont. 

2131.3 

341.6 

1492.9 

3.04 

7069.5 

1359.9 

cont. 

2068.7 

336.0 

1436.6 

2.88 

7173.7 

1348.7 

cont. 

1993.2 

330.6 

1366.9 

2.73 

7300.2 

1338.0 

cont. 

1916.3 

326.0 

1295,2 

2.59 

7430.1 

1327.7 

cont. 

1810.0 

321.8 

1223.4 

2.46 

7560. 2 

1317.6 

cont. 

1766. 0 

318.2 

1153.5 

2.34 

7687.3 

1307.9 

cont. 

1695.5 

314.9 

1086.5 

2.23 

7609.9 



Tempo! 11(88) Rendezvous Titan m E/Centaur/Sert III 3-Year Missions 


AD 

r mln 

r 

max 

r 

com 

9 

com 

*mtn 

0 

max 

AX 

6 

1 

LD 


S 

m 

o 

m 

P 

m . 
net 

P t 

V 

m 

0 

1.010 

3.168 

.96 

89.4 

70.1 

131.3 

363.4 

-15.53 

28.50 

7/29.7/85 

1144.8 

cont. . 

2272.7 

428.3 

1539.0 

5.45 

6836.4 

-10 

1.007 

3.191* 

.92 

91.9 

67.1 

134.0 

355.1 

-15.07 

28.50 

7/30.7/85 

1133.8 

cent. 

2172.3 

419.2 

1448.6 

5.43 

7001.6 

-20 

1.004 

3.212 

.89 

94.9 

63.1 

136.7 

347.4 

-14.19 

28.50 

7/31.2/85 

1123.3 

cont. 

2063.1 

410.6 

1348.8 

5.35 

7183.1 

-30 

1.000 

3.233 

.86 

98.6 

58.4 

139.3 

340.3 

-12.92 

28.50 

7/31.2/85 

1113.3 

cont. 

1052.0 

402.8 

1246.3 

5.23 

7369.7 

-40 

.994 

3.253 

.83 

103.0 

53.4 

141.8 

333. 9 

-11.36 

28.50 

7/30.7/85 

1103.8 

cont. 

1843.0 

395.8 

1145.1 

5.08 

7554.9 

-60 

. 930 

3.292 

.79 

114.3 

48.3 

145.7 

324.2 

-7.68 

28.50 

7/27.4/85 

1087.1 

10C9.1 

1626. 9 

372.7 

954.3 

4.69 

7030.2 

-70 

.971 

3.311 

.78 

121.2 

46.1 

147.2 

320.6 

-5.89 

28.50 

7/24.9/85 

1079.6 

1051.2 

1527.9 

3G2.0 

867. 1 

4.47 

8105.9 

-80 

.901 

3.331 

.77 

129.0 

43.3 

148.3 

317.8 

-4.13 

28.50 

7/21.8/85 

1072.7 

1035.2 

1437.3 

353.6 

785.7 

4.25 

8269.4 

-90 

.919 

3.352 

.78 

137.3 

40.2 

148.9 

316.2 

-2.33 

28.50 

7/17.0/85 

1066. 5 

1021.0 

1355.3 

347.5 

710.5 

4.02 

8419.5 

-100 

.934 

3.373 

.79 

145.5 

37.0 

148.8 

316.0 

-0.30 

28.50 

7/12.9/85 

1061.6 

1000.4 

1282.5 

344.2 

641.2 

3.81 

8554.6 

-110 

.915 

3.395 

.83 

152.4 

34.1 

147.5 

318.4 

2.38 

28.50 

7/5. 5/85 

1059. C 

1003.1 

1221.4 

346.2 

578.0 

3.60 

8669.6 


o 

oo 


2-Year Missions 


0 

.825 

2. %35 

.95 

89.4 

54.0 

128.8 

414.2 

7.14 

28.50 

6/6.7/86 

832.8 

cont. 

1346.8 

462.8 

575.2 

^.45 

8435.X 

-10 

.792 

2.465 

.92 

91.9 

57.5 

128.7 

416.6 

9.21 

28.50 

5/27.6/86 

832.9 

cont. 

1297.9 

459.7 

519.6 

5.43 

8344.6 

-20 

.761 

2.497 

.89 

94.9 

58.9 

129.9 

417.5 

10.27 

28.50 

5/1 D. 1/86 

831.4 

cont. 

1293.9 

454.7 

466.2 

5.35 

8655.3 

-30 

.731 

2.531 

.86 

99.6 

59.0 

169.0 

417.6 

10.75 

28.50 

5/11.6/96 

828.9 

cont. 

1171.0 

449.5 

415.1 

5.23 

8^65.4 

-40 

.701 

2.566 

.83 

103.0 

59.3 

169.3 

417.3 

10.89 

28.50 

5/4 . 8/9G 

825.7 

cont. 

1114.7 

441.5 

366.5 

5.09 

8873.6 

-50 

.671 

2.602 

.81 

109.2 

57.3 

167.7 

416.8 

10.85 

28.50 

4/28.5/86 

822.0 

cont. 

1060.6 

434.0 

320.5 

4.99 

8979.0 

-60 

.612 

2.639 

.79 

114.3 

56.0 

167.3 

416.2 

10.73 

28.50 

4/22.6/86 

817.9 

cont. 

1009.9 

426.2 

277.5 

4.69 

90*0. 7 

-70 

.612 

2.676 

.79 

121.2 

54.7 

166.9 

415.6 

10.58 

28.50 

4/17.2/86 

813.3 

cont. 

959. 9 

413.7 

237.3 

4.47 

9179.4 

-80 

.583 

2.714 

.77 

129.0 

53.4 

166.6 

415.0 

10.43 

28.50 

4/12.1/96 

809.4 

cont. 

913.6 

410.0 

200.0 

4.25 

9271.6 

-90 

.555 

2 751 

.79 

137.3 

52.1 

166.3 

414.3 

10.32 

29.50 

4/7.:i/"fi 

803.2 

cont. 

870.1 

401.8 

165.5 

4.02 

9360.3 

-100 

.527 

2.788 

.79 

145.5 

50.9 

166.1 

413.7 

10.26 

28.50 

4/2.9/90 

797.6 

cont. 

829.4 

393.7 

133.6 

3.81 

9444.2 




Terapel 11(88) Reodacvous ShuttWTraiiatago/Sert III 


AD 

r mtn 

r 

max 

r 

com 

0 

com 

*mln 

# 

max 

AA 

6 

1 

LD 

0 

1.018 

3.768 

.98 

89.4 

50.2 

104.2 

400.3 

-3.51 

28.50 

6/20.7/84 

-to 

1.018 

3.777 

.89 

94.9 

50.1 

109.9 

383.9 

-2. 88 

28.50 

6/22.7/84 

-40 

1.018 

3.789 

.83 

103.0 

50.0 

115.9 

368.0 

-2.32 

28.50 

6/24.7/84 

-60 

1.016 

3.803 

.79 

114.3 

50.0 

122.1 

353.2 

-1.82 

28.50 

6/26.7/84 

-80 

1.016 

3.819 

.77 

129.0 

49.9 

128.4 

339.7 

-1.37 

28.50 

6/28.7/84 

-100 

1.016 

3.836 

.79 

145.5 

49.9 

134.5 

327.5 

-.95 

28.50 

6/30.6/84 

-120 

1.016 

3. 855 

.88 

156.1 

49.9 

139.9 

316,6 

-.52 

28.50 

7/2. 5/84 

-140 

1.016 

3.877 

1.04 

149.8 

49.8 

143.8 

306.9 

-.04 

28.50 

7/4.3/84 

-180 

1.016 

3.902 

1.27 

134.9 

49.7 

146.9 

298.1 

.57 

28.50 

7/6.1/84 

-180 

1.016 

3.932 

1.56 

119.4 

49.4 

150.2 

290.0 

1 42 

28.50 

7/8.0/84 

-200 

1.015 

3.975 

1.90 

104.7 

48.4 

152.9 

282.2 

2.76 

28.50 

7/10.2/84 

-220 

1.015 

4.032 

2.27 

90.6 

45. C 

149.3 

273.2 

4.20 

28.50 

7/14.3/84 

-240 

1.014 

4.084 

2.64 

77.1 

41.2 

146.9 

264.5 

5.30 

28.50 

7/18.7/84 

-280 

1.014 

4.126 

3.01 

63.9 

35.9 

146.4 

256.7 

5.64 

29.50 

7/22.6/84 

-280 

1.012 

4.159 

3.35 

51.0 

30.6 

147.2 

249.9 

5.47 

28.50 

7/25.6/84 

-300 

1.011 

4.186 

3.66 

38.4 

26.1 

148.7 

244.3 

5.19 

28.50 

7/27.8/84 

-320 

1.009 

4.211 

3.92 

26.1 

22.3 

150.4 

239.8 

5.01 

28.50 

7/29.2/84 
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Tempel 11(83) Rendezvous 
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APPENDIX B 


The equations for the partial derivatives av /5v o and 9v /S6 are 
derived from Equations (5) and (6) of Section IV and are listed below: 

2 2 2 2 
9v v r 9 w(3+2pw-w )(9 p/9 V ffl ) + (3p + 2p w -3pw -4w)(9\v/9v oo ) 

— £- = — — < p ' _ — 

dv v l 9 v / 1 1 v 

® 6 2.y(l+pw-w )(2+pw) 

2 2 2 

9v v (3p + 2p\v-3pw -4\v) ^ 

g _ o v ^ 9w 

96 “ , 2 S6 ’ 

2v /(l+pw-w )(2 + pw) 

S 


where 


9p/9v = l/v , 

CD O 


9w 

_ 1 

f 1 . a .v , r 

7, x - - 

)^2L + (£ + y ) 

dv 

00 

2 

i — — + — *- -f 

t2 9v 9v L 

00 00 

U+ / 9 o 

A 2 7xV4 + S“ 

) 9v v2 y / 

00 


+ )]/(2w-p/2-y)j , 

CO 

T TSy + x(9x/9 6 ) +4s(9 s/9 6) + fB + y \ ill/ (2w -p/2 - y)| 

96 2 19 6 L 96 0 / 277-j >56j 1 


9 
du 


2, Zx/4 + s‘ 

9 s/9 6 = cos (6 - i), 

a y/ s v oe = [(p/2) 9p/9v ao -9x/9v e6 ]/2y, 

9y/9 6 = - (9x/d6)/2y, 

.2 


iL = I T (q/2)(9q/5u) + (p/3) 2 (Sp/9u) + |SL *| [ / (q/2) 2 + (p/3) 3 + q/2 1 
6L ./(q/2?^ (p/S?” ^ J 


-2/3 


,ip/2)C9q/9u) + ( P /3) 2 ( a p /9„) _ |* J _ q/JJ 7 
/( q/2f + (p/3) 3 


-,—2/3 


with u = v or 6, 
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Sq/9v oo = 2ps 2 (l -s 2 )(9p/3v oo ), 

3q/36 = 2p 2 s(l -2s 2 )(3s/3 6), 

2 

dp/dv oo =2ps (Sp/dv^), 
dp/3 6 = 2s (p 2 +4) (3 s/3 6). 
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